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Abstract 


A new deterministic mathematical model for the transmission dynamics of Middle 
East Respiratory Syndrome Coronavirus (MERS-CoV) is proposed and fully analyzed. 
The presented model exhibits a unique endemic equilibrium and there is no infection free 
equilibrium due to constant influx of latent immigrants. An invasion threshold parameter is 
derived and stability analysis of the full model and its two special cases is carried out. The 
impact of quarantine and isolation is assessed via threshold analysis approach, while the 
impact of immigration on the disease prevalence is discussed. Indeed, we showed that MERS- 
CoV can be controlled by quick isolation or monitoring close contacts and quarantining of 
suspected latent immigrants. Further, numerical simulations of the model reveal that the 
disease can be contained if these preventive measures are combined with high reduction of 
immigration rate. 


Subject Classification: (2010) 34D23, 93C15, 93D20 
Keywords: Middle East Respiratory Syndrome, Isolation, Quarantine. 


1. Introduction 


Middle East Respiratory Syndrome Coronavirus (MERS-CoV) belongs 
to the subfamily Coronavirinae in lineage C of the genus Betacoronavirus. 
MERS-CoV is one of the four Coronavirus groups of which Severe Acute 
Respiratory Syndrome (SARS-Cov) belongs [9]. In September 2012, The 
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World Health Organization (WHO), reported the First cases of MERS- 
CoV [26], a novel respiratory disease, initially localized to Middle East 
countries. It has a high potential for transmission to close contacts, as 
observed in communities with sporadic exposure, health care facilities, as 
well as families having contact with infected members [10]. Coronaviruses 
are wide spread in Bats around the world but can also be found in 
many other animal species as well [20]. The first index was reported in 
Kingdom of Saudi Arabia (KSA) on 13th June, 2012, when unknown 
coronavirus isolated from a 60yr old man admitted in hospital with a 7 - 
day history of acute pneumonia, fever, cough, expectoration, short breath 
and subsequent renal failure which later resulted in fatal outcome [28]. 
Shortly thereafter, on 14" September, 2012, the United Kingdom Health 
Protection Agency (HPA) Imported Fever Service was notified of a case 
of unexplained virus, later named human coronavirus England 1. It has 
been isolated from another patient with history of severe respiratory 
illness in an intensive care unit (ICU) of London hospital. The index has 
been transferred from Qatar and had a history of travel to Saudi Arabia 
[3]. The origin of the disease was traced back to an earlier time - April, 
2012, when health care workers in an ICU of a hospital in Zarqa, Jordan 
were confirmed to have been contracted with the novel coronavirus [9]. 
From June, 2012 to 25" September, 2013, WHO reported total of 133 cases 
with 55 deaths from 9 countries [20]. Up to 14" March, 2014, there were 
477 globally reported cases and alerts including additional 7 countries 
[15]. Globally, from September 2012 to 17" April, 2013, WHO has been 
informed of a total of 243 laboratory-confirmed cases of infection with 
MERS-CoV, including 93 deaths [25]. It has been observed that all cases 
outside Middle East are traced back to the region either by travel to or 
from the affected areas or by direct or indirect contact with others who 
have a travel history to Middle East [1]. 

The epidemiological data available suggest that the infection is 
primarily zoonotic in nature, with possible cross-species transmission 
between humans and animals (Bats, Camels). Some of the reasons advanced 
are: Close phylogenetic similarity between viruses isolated from humans, 
bats and camels; presence of viral neutralizing antibodies in dromedary 
camels for the past 2 decades from (EGP, JO, KSA, OM, Canary Island), 
and none in other livestocks [1, 13, 16, 21, 22]; close association between 
camels and humans as source of food, sporting and as pets. However, 
most of the reported cases have no close contact with camels or bats . 
Therefore, the direction of transmission from humans to animals or vice 
varsa, host, reservoir e.t.c. remain unknown [1, 16, 17]. 
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The established clinical symptoms of MERS-CoV include fever, 
cough, shortness of breath, acute pneumonia, expectoration, hemoptysis 
often followed by renal failure [2, 20, 28]. Gastrointestinal symptoms like 
vomiting and diarrhea are also common [2]. The incubation period was 
estimated to be 5.2-5.5 days and the serial interval (the time between the 
successive onset of symptoms in a chain of transmission) was 7.6 days 
[2]. According to Cauchemez et al, [7], the estimated basic reproduction 
number obtained from all confirmed and probable of human cluster cases 
from Middle East (up to Aug 8, 2013) and the first 133 cases [20], were in 
the range 0.81.3 and 0.6 respectively. The implication of the former result 
is that self-sustaining transmission in the absence of control measures 
may occur. The average age of the first 133 reported MERS-CoV cases was 
52 years with male having higher case-fatality risk of 52% compared to 
women 24% [20]. 

Clinically, whenever there is an outbreak or re-emergence of new or 
an old infection in a community/country, in the absence of vaccine or 
treatment, isolation and quarantine of individuals with clinical symptoms 
or suspected to be exposed to the pathogens have proved to be effective 
ways of containing the disease [11, 13]. The World Health Organization 
is concerned, among others, with the gaps in understanding the degree 
of transmissibility between people, possibility of “super spreaders” and 
potential for sustainable person-to-person transmission and spread [10]. 
However, several mathematical epidemic models incorporating one of 
such control measures or both exist in the literature, see for instance, [11, 
18, 23, 24, 27]. For the 2015 epidemic outbreak of MERS-Cov in the Republic 
of Korea, two deterministic models for the transmission dynamics of 
MERS-Cov with and without control measures were developed in [27]. 
The basic reproduction number, R, was estimated to reach up to 4.422. 
The numerical analysis reveals that lack of self-protection sense and 
targeted control measures were the reasons of the outbreak spread quickly. 
However, it was reported that strengthening self-protection ability of 
susceptible and quickly isolating or monitoring close contacts are effective 
measures to control the disease. Furthermore, partial correlation analysis 
shows that the infectivity and proportion of the asymptomatic infected 
cases have much influence on the disease spread. In the same vain, another 
deterministic model was designed to analyze the MERS-CoV outbreak in 
the Republic of Korea [18]. The proposed model explicitly incorporates 
superspreading events and time-dependent transmission and isolation 
rates. The superspreaders are those who transmit the virus to more than 
20 patients and have underlying respiratory diseases with a severe cough. 
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The impact of the timing of control measures associated with a reduction 
of the transmission rate and diagnostic delays on the outbreak size and 
duration was assessed. It was observed that early interventions within 1 
week after the epidemic onset is a promising means to reduce the size and 
duration of the MERS-CoV epidemic. 

In a similar note, an SEIQJR - SARS model with both quarantine 
and isolation strategies was presented in [11]. The analysis of this model 
reveals that a reduction in the contact rate between susceptible and infected 
individuals by isolating the latter is a critically important strategy that 
can control SARS outbreaks with or without quarantine. It was reported 
that an optimal isolation is more effective than sub-optimal isolation 
and quarantine together. And that an optimal isolation combined with a 
highly effective screening programme at the points of entry would lead to 
a community-wide eradication of SARS. 

It is well known that migration and immigration greatly increase 
the spread of many infectious diseases at a regional, national and global 
scale [4, 12]. Simple models for disease transmission incorporating 
immigration of infective individuals was presented in [4]. The models 
exhibit a single endemic equilibrium that is asymptotically stable. It 
was revealed that there is no disease free equilibrium in the presence of 
immigration of infectives. For HIV transmission model a considerable 
reduction of infectives was suggested by screening and quarantining of 
infectives. While both immigration/migration terms were incorporated 
into all sub-population compartments, susceptible and infected, of two 
types of well-known heterogeneous epidemic models: multistage models 
and multi-group models for HIV/AIDS and other sexually transmitted 
diseases (STDs). It was shown that, the disease always becomes endemic 
in the population and tends to a unique asymptotically stable endemic 
equilibrium when migration or immigration into infected sub-population 
is present [12]. 

In this paper, a disease transmission model for MERS-CoV epidemic 
incorporating immigration of susceptible and latently infected individuals 
is proposed to assess the impacts of quarantine and isolation strategies in 
controlling the disease. As in the previous studies, the main result is that 
our model exhibit a unique asymptotically stable endemic equilibrium 
and there is no disease free equilibrium when immigration into latently 
infected sub-population is present. While such equilibria exist for the case 
when either all the immigrants are susceptible or the immigration rate is 
negligible. The impact of quarantine and isolation is assessed in relation 
to the relative infectiousness of quarantined and isolated individuals. 
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Further, eradication of MERS-CoV is feasible by reduction of immigration 
rate of latently infected individuals. The proposed model extends in some 
sense the aforementioned models in the literature as follows: 


(i) In the presented model, demographic features are incorporated 
and vital mathematical analysis is carried out, while these were not 
provided in MERS-CoV model proposed in [27]. 


(ii) In [18] there were only two sources of infection, here we consider 
four sources of infection. 


(iii) Although, immigration into both susceptible and _ infective 
compartments is common to [4, 12] and our model, in the latter 
we used saturated incidence formulation when the former used a 
simple mass action formulation. 


(iv) For the SARS model in [11], quarantined individuals are 
asymptomatic infectives who develop symptoms and then move to 
the isolated class. In our model, some quarantined individuals do 
not show symptoms and become susceptibles while others recovered 
by either treatment or acquiring natural immunity and move to 
recovered compartment. 


The remaining part of the paper is organized as follows. In section 
2, we describe the mathematical model. In section 3, rigorous analysis 
of the model is provided including its basic properties, existence and 
stability of the model equilibria, and the impact of quarantine, isolation 
and immigration rate are discussed. Numerical simulations supporting 
the theoretical results are provided. In section 4, we provide a concluding 
remarks. 


2. Model formulation 


The susceptible human population is generated via recruitment by 
birth (at a constant rate II) and immigration (at a rate aq); where 0 <q < 
1 account for the fraction of the inflow of migrants into the community 
(by air or road) who are susceptible, and by recovery of quarantined 
individuals at a rate o(1 — p); where 0 < p< 1 isa fraction of quarantined 
individuals tested with no clinical symptoms of the disease. This 
population is decreased following infection with MERS-CoV, which can 
be acquired via effective contact with infectious human (at a rate A), where 
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_ AQ,L+l+n,Qt+n,J) 
N 


A 


Here the parameter f is the effective contact rate (contact capable 
of leading MERS-CoV infection). Furthermore, 0 < n, < 1,0< 7, <1 and 
0 < n, < 1 are modification parameters accounting for the assumed 
reduction in infectiousness of individuals in the latent (L), quarantine (Q) 
and isolated (J) classes, in comparison to infectious individuals in I class. 
Natural death is assumed to occur in all human compartments at a rate LU. 
Thus, the rate of change of the susceptible population is given by 


= = Il+agt+o(1-p)Q-(A+ ws. 


The population of individuals in the latent class is increased by 
immigration at the rate o(1 - q) and by infection (at the rate A), and is 
decreased by the development of clinical symptoms (at a rate 1), 
quarantine (at a rate t,) and natural death, so that 


a(1—q)+AS—(t, +7, +m)L. 
dt 


dL = 


The population of infectious individuals with clinical symptoms 
of MERS-CoV in I class increases following the development of clinical 
symptoms by individuals in latent class (at the rate t,). This population 
is decreased by isolation (at a rate y), recovery (at a rate 6,), MERS-CoV 
induced mortality (at a rate v,) and natural death, this gives 


dl 
Ti 7,L-(yt+v,+6,+w)I. 

The population of quarantine individuals increases following the 
quarantine of individuals in the latent class (at the rate t,). This population 
is decreased by recovery (at a rate 0) and natural death, so that 

dQ 
— = 1,L-(o+ ; 
cf gh-(a+4)Q 

The population of individuals that are isolated (J(#)) is generated by 
the isolations of infectious individuals with clinical symptoms of MERS- 
CoV (at the rate 7). It is decreased by recovery (ata rate 6,), disease induced 
death (at a rate v,) and natural death. Hence 
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(utv2)J 
opQ 


Figure 1 
The flow diagram of model (1). The parameter, 4 = a 


d 
“ = yI-(v,+6,+ yw). 

Finally, the recovered population is generated by the recovery of 
individuals in I, J] and Q classes at the rates 6,, 6, and op, respectively, and 
decrease due to natural death. This gives 


dR 
ae 0,1 + 6,1+opQ- u,R. 


In summary, the MERS-CoV transmission model is given by the 
following system of non-linear differential equations (a flow diagram of 
the model is given in Figure 1 and the associated variables and parameters 
are described in Table 1, respectively). 

dS AL +1+7,Q+nJ) 

a Il+aq+o(1—p)Q N 
De. wieaja AL +1+7,Q +75) 

dt N 
dl 
dt 
dQ 
dt 


S-uS, 


S Tot + OL, 
= 1,L-(y+v,+6,+m)I, (1) 


= 1,L-(o+ WQ, 
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d 

a = yI-(v,+6,+)J, 
dR 

rs = 6,1+6,J+opQ-uR, 


with nonnegative initial conditions 5S(0)2=0,L(0)=0,1(0) =0,Q(0) = 0, 
J(0) = 0,R(0) = 0. 

Since the last compartment (R) does not appear in the first five 
equations of model (1), it is convenient to discard the last equation of 
system (1) for theoretical analysis. 

The equation for the total population which is obtained by adding all 
the equations of (1) is given by 

dN(t) 


<= T+a-aN()-v,10)-v JO: (2) 


3. Analysis of model (1) 
3.1 Basic properties 
The basic properties of the model can now be investigated. 
Lemma 1: The solution (S(t),L(£),I(t),Q(t),J(£)),R(£)) of model (1) with positive 
initial conditions, exists for all t > 0 and is unique. Furthermore, S(t)>0, 


L(t) > 0,1(t) > 0,Q(f) > 0 and J(t) > 0 for all t = 0. 
Using Lemma 1 on equation (2), we have the following result. 


Lemma 2: The biologically-feasible region of model (1) is 
T= {(S¢),L(t),1(t), Q(), I(t), R()) € RE : S(t) + L(t) 


+1(t)+ Q(t) + J() + R(t) < U 2) 
u 


which is positively-invariant and attracting. 


Proof : Adding the equations of the model (1) gives 
dN(t) 


— = +a-yN(t)-v,I¢)-v,J@), 


so that 
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dN(t) 


OA = H+a-pN(). 
ri a — uN (t) 
It follows that ee <0 if N(t)> ee, Thus, a standard comparison 
theorem can be used to show that 
Nit) < Ne + tee], 


Il+a@ 


In particular, N(t)< if N(O)< —. Thus Tis positively 


Il+a@ 


invariant. Furthermore, if N(t) > , then either the solution enters T 


in finite time or N(t) approaches !!*@ and the infected variables L, I, Q 
u 


and J approaches zero. Hence, IT is attracting. That is, all solutions in R; 
eventually enter Tl. 


Table 1 
Description of parameters of model (1) 


Parameter Interpretation 
II Recruitment rate of susceptible humans by birth 
q Fraction of immigrant humans 
a Recruitment rate of humans by immigration 
Oo Removal rate from quarantined 
B Transmission coefficient 
Ny Ny Ns Modification parameters for reduction in infectiousness of 
latent, quarantined & isolated individuals, respectively 
Lu Natural death rate in all classes 
Tp Ty Progression rates to infectious / quarantined classes, 
respectively 
Y Isolation rate 
Vi, V> Disease-induced death rates for infectious and isolated 
individuals , respectively 
6, 5, Recovery rates for infectious and isolated individuals , 
respectively 


p Fraction of quarantined humans 
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3.2 Existence and stability of equilibria 
It should be noted that model (1) has no disease free equilibrium (DFE) 


A AA RA 


be the equilibrium solution of system (1). Then setting the right hand side 
of (1) to zero we obtain 


a K,K,(11 + @q)+ aor, (1—p)(1—4) 
(1+ A)[r + (t+ (0+ w))+7,0(U+ Ap)’ 
K,[A(1+@)+ wa(1—-q)] 


[= —, 
(ut ANT, u+(7,+ WY(o+ w)]+7,0(u+ Ap) 
ve 7,K,[A(1+ @)+ ya(1-q)] (3) 


K (ut Aru t(r, +W(o+ W]+7,0(u+ Ap)} 
7, [A(I1+ @) + wa(1—q)] 
(ut A)[r,u+(r, + u(o+ W]+7,0(u+ Ap)” 
- T,7K, [A(I1 + a) + uwa(1—q)] 
KK (ut Ac,ut(c, + (ot w+ 7,0(u+ Ap)} 


> 


where K, =7,+7,+u,K, =y+v,+6,+u4,K,=o+y,K,=v,+6,+y, 
R - BI +n,L+n,Q+n,J) 
N 


(4) 


and 
N= StLel +04}, (5) 
Substituting (3) and (5) into equation (4) yields after algebraic 
manipulations (6). 
A{K,K, (11 +g) + aor, (1- p)(1-4)} 
(u+A)[t,u+(z, + w)(ot+ w)]+7,0(u + Ap) 6) 
4 7,K,4-£) 7,7K,(A-7,£) =0 
K KR 


2 


+ ®K,(A-7,f) 


+7,(A 7,8) + 


with 
— di+a)+ne-g) 
(ut+A)[z,u+(7,+ (ot W)]+7,0(u+ Ap) 


MERS-COV WITH LATENT IMMIGRANTS 11 


Thus, the endemic equilibria of system (1) correspond to the positive 
values of (3) obtained by using the positive solutions of (6). 
We define an invasion threshold by 
_ BUI + a@)[t,K,(K, +7,7)+K,K,(7,K, +7,7,)] 


R, = ose eee (7) 
K,K,K,K, (11+ aq) 


Thus, we obtain the following results. 


Theorem 1 : Model (1) has: 
(i) a unique endemic equilibrium for all values of R, if q € [0,1) and a> 0, 
(ii) a unique endemic equilibrium if R,> 1 and either q=1 or a=0, 
(iii) no endemic equilibrium if R, <1 and either q = 1 or a= 0. 
Proof : Re-arranging and simplifying equation (6) gives the following 
quadratic equation for 4 
al? +bA+c = 0, (8) 
where 
a = (I1+qa@)[K,K,(K, +7,)+7,K,(K,+y)], 
4 
b = [Ka +aq\1—-R,)+a(1—q)ut{t,K,[K,o(1—p)+K,u)] 
i=1 
4K, K jt, +K,)4+7,Kcu)|, (9) 


c -apu (1 ve g@)[K,K,(7,K, + 1,7) + 7,K,(K, + 73Y)1, 


The coefficient a of equation (8) is always positive and c is always 
negative (noting that ~@>0 and q é€ (0, 1)), whereas b is positive if R,< 1. 
If R, > 1, b can also be positive or negative. Thus, model (1) has a unique 
endemic equilibrium for any sign of b by Descartes rule of sign. This 
proves (i). Next, if q = 1 or a = 0, then the coefficient ** and equation (8) 
reduces to 


aA+b =0, (10) 


with either 


a 


C+ e)( 5 KAR tooo KR Fy), (11) 


b = [IK (1+ a)-R), 


or 
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a 
Il 


NK KAK +) +7, K+ 7), (12) 


4 
[[K,d-R,), 
j=l 


> 
ll 


respectively. It is clear that, in each case the solution A= -" of (10) is 
positive if R, > 1 and negative when FR, < 1. Moreover, for Rk, = 1, the 
coefficient b of (10) is zero and so, the coefficient a of this equation must 
be zero since 2 cannot be zero. Hence, equation (10) has no solution 
whenever FR, = 1. Consequently these prove (ii) and (iii). oO 


It follows from (7) that if either q = 1 or a = 0, then the invasion 
threshold, R, reduces to 


pe = Blt Ks(Ky +7) + KKK +7) a3) 
K,K,K,K 


1°27 "3° "4 


The uniqueness of the endemic equilibrium, € allows us to 
numerically show its global stability as depicted in Figure 2, since the 
analytical proof is quite involved. 


3.3 Immigration cases 


As mentioned earlier, when all the immigrants are susceptible, 
q = 1 or the immigration rate is negligible (@ = 0) the threshold parameter 
is as defined in (13). Furthermore, for each of such scenarios a disease free 
equilibrium exists. 

For the first case, in the absence of infection, model (1) with g = 1 
exhibits a locally asymptotically stable disease free equilibrium given by 


Il+a 


= CLT -[ 0,0,0,0} 


In fact, the Jacobian matrix evaluated at &, denoted by a , is 


= =f, B o(l-p) —£n, 


O° fy Hk, =B- fay. . Bi, 

3, =| 0 Ci -K, 0 Oo |. (14) 
0 tl, 0 -K, 0 
0 0 y 0 -K, 
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The eigenvalues of this matrix are: —u, and the eigenvalues of the 
matrix J g Which is obtained by deleting the first/last rows and the first/ 
last columns of J,- Then, the characteristic polynomial of this matrix is 


At +A? +B? +CA+D=0, (15) 

with 

A = K,+K,+K,+K,- 7, 

B = K,K,+(K,+K,)(K, +K, — Bn,)+K (K, — £n,)- Bt, +7,7,) 

K K 

Br{K, ID) KK, kof BOK, st) 

KR A, K,K, 

+K,K,K, — B(7,K,K, + 7,K,) 
= K.K,K,(1-R,,)+K,K,(K, +K,)Q—-R,,)+K,K,K, 


=P KK + 7,K,) 
D= K,K,K,K,(1-R;). 


C= KK,K, fs 


The Routh-Hurwitz criteria [19] for the local asymptotic stability in 
this case are A>0,C>0,D>0 and ABC-—C*-A’D>0. The first three 
conditions are satisfied if R; <1 which implies that R, <1 and R, <1 
so that Br, <K,K, and £7,<K,, respectively. The last condition holds 
since all terms remaining in the expansion after cancelation are positive. 
Thus the disease free equilibrium, €j of model (1) with q = 1 is locally 
asymptotically stable when R, <1. 

It is instructive to note that if R; =1, the coefficient D of the 
characteristic equation (15) is zero and so, J, has a zero eigenvalue. This 
indicates the occurrence of a bifurcation when R; =1 which will be 
investigated using Centre Manifold theory as described in Castillo-Chavez 
and Song [6, Theorem 4.1] in Subsection 3.4. 

The epidemiological interpretation of this result is that MERS-CoV 
will be eliminated if the initial sizes of the infected sub-populations of 
the model are in the basin of attraction of the DFE (€!). However, disease 
eradication is independent of the initial sizes if the DFE is globally 
asymptotically stable. The global stability result is established in the 
following theorem. 


Theorem 2: The DFE (€) of model (1) with q = 1 is globally-asymptotically 
stable (GAS) when R, < 1. 
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Proof : The proof is based on the approach in [1]. Using the notation in [5], 
for system (1) with q=1,X =S,Z =(L,I,Q,J), and X' = (7,0). Then, 
the equation of uninfected is - 

dx 


dt 


= F(X,Z)=Il+a+o(1-p)Q-AS—uS, 


and the infected subsystem is 
ASC. T YL 
aZ G(X,Z) = 7,L-(y+v,+6,+ 1 
dt t,L—-(o+ “)Q 
yl-(v,+6,+u)J 
When L=I=Q=]J=0 (ie., Z =0), then F(X,0)=(1+a@)—uS, and 
Bbn,-K, B Bn, Pn, 


-K, 0 0 
A= Det oj=|" 2 
: 0 =k. 0 


0 y 0 —-K 


4 


is obviously an M-matrix (a matrix whose off diagonal elements are 
nonnegative). 

Next, we have 

G(X,Z) = AZ-G(X,Z), which implies that 


ane{r-S)ei(1-F} rn o[1-S)eni(2-S)} 
G(X,Z)= 0 


0 
0 


It is clear that CCS Z) =0, since 0 < S < N. It is also clear that 
* T+ ; F i ae aX _ 
X= [72,0] is a global asymptotic stable equilibrium of [= F(X,0). 


Thus, by using a result in [5], we conclude that € is globally-asymptotically 
stable when R; < 1. im 


It follows from Theorem 1 that, model (1) with q = 1 has a unique 
endemic equilibrium if R; > 1. 
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Now, let €¥ and a be the disease free and endemic equilibria of 
model (1), respectively with negligible immigration rate (@ = 0), then we 
obtain the following results. 


Theorem 3 : The unique DFE, €¥ of model (1) with a = 0 is globally- 
asymptotically stable if R; < 1, otherwise unstable. 


Theorem 4: The unique endemic equilibrium, E of model (1) with a = 0 is 
globally-asymptotically stable if R; > 1. 

Hence, the summary of the stability results of the model equilibria is 
given in Table 2. 


3.4 Bifurcation analysis 


Here we investigate the existence of bifurcation of model (1) using 
Centre Manifold theory since the Jacobian matrix J, for the case when 
all the immigrants are susceptible has a simple zero eigenvalue when 
R; = 1 with other eigenvalues having negative real parts as presented in 
Subsection 3.3. In order to apply such a theory, we make the following 
change of variables S=x,,L=x,,1=x,,Q=%,,] =X,, so that we use the 
vector notation x = ie a Fe Pree )’. Then system (1) with q = 1 can be 


d. 
written in the form <= F(x) with f =(f,f-fy- Sirf.) » such that 


X, +X, +7,X,+7,X 
k, = Il+a+o(1-p)x, BUX, + %5 ¥Th%s + Ms a) ae 
X, +X, +X, +X, +X, +X, 


‘ XL, +X, +9,X, +7,X 
i, = AUR ra al dy (1, +7, +1)X,, (16) 
KX, +X, +N, +X, +X +X, 


= Fx, ~ (FEV FO, FB) 
= 1,X,-(O+p)Xx,, 
ty = ¥X,—(V, +0,+ WX: 

The Jacobian matrix of system (16) at €, is the same with J, as presented 
in (14). If Bis taken as a bifurcation parameter, we obtain from R; = 1 that 
K,.K,K,K 


B = f= 6 lie sis: Nader ; 
7,K,(K, +7,7)+ K,K,(7,K, +7,7,) 


It can be verified that the Jacobian of (16) at B = f has a right 
eigenvector associated with zeroeigenvaluegivenby u =(u,,U,,U,,U,,Us) , 
with 
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u, = ely, and u, => (17) 


Similarly, the left eigenvalue of J,(€, 0) associated with zero eigenvalue 
at B = Bis given by 0= (0,040,040, )’, where 


B Bn Bn 
v, =0,0, =v, > 0,0, “KK -K,)v,,0, = Z 2”,,0, =— 


3 4 


(18) 


To establish the existence of bifurcation, we need to compute a and b 
as follows. Let F, be the k” component of F, then as defined in [6] 


oF, | 
Eras é) 


k,i,j=1 


ues 


k,i,j=1 xX 


a 


(19) 


> 
Il 


Now, for system (16) the associated non zero partial derivatives of F 
at the disease free equilibrium are as follows: 


OF, Bonu OF, fun, +1) 


x; +a’ Ox,X, Il+a (20) 
OF _ Buatm) &F _ Bunn) 
Ox,X, Mea” Ox,X, +a 


It follows from (19) and (20) that 
7 Belts fan, AnD) om tm) + 1) 


(21) 
+a K K K 


2 3 4 


which is negative, and for the sign of b, it is associated with the following 
non zero partial derivatives of F 


aE OE OE OE 
mop" &x,o8 oxo Oxop 7 
2 3 4 5 
OF OE OF OF 
— Ui 2 =1, ; =o, ; =). 
ax,08 ax,08 ° ox,08 ” Ox,0p 
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Table 2 
Summary of the stability results of the full model and its two special cases. 


Equilibrium Invasion threshold Stability 
e R,S1 GAS 
& 21 LAS/GAS 
ce R;<1 LAS/GAS 
é, Ri>1 GAS 
E Ris GAS 


It follows from the expression of b in (19) and the derivatives in (22) 
that 


Tt i 
be [ne ae >0. (23) 


2 3 4 


We obtain from the computations of a and b above that a is negative 
and b is positive. Thus, using Theorem 4.1 (iv) in [6], we establish the 
following result 


Theorem 5 : The system of equations (1) with q = 1 exhibits a transcritical 
K,KyK3Kg 
7 K3(Ka+137)+K Ka (mK3 +1972) aid 


bifurcation at €, if the bifurcation parameter B= 
equivalently R; = 1. 


3.5 The impact of quarantine/isolation/immigration 


In this section, the impact of quarantine and isolation on the 
transmission dynamics of MERS-CoV is investigated via a threshold 
analysis approach on the invasion threshold parameter, R,. This is carried 
out in terms of the parameters associated with quarantine of latent 
immigrants 7, and the isolation of individuals with disease symptoms y. 
More precisely, the threshold R, is considered as a function of t, and y. 
For the case of the quarantine of latently infected individuals, the rate of 
change of , with respect to 7, is as follows: 
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OR, Bn, #IK,K,(7,K, +7,7,)+7,K,(K, +7,7)] 


- 5 ‘ (24) 
Ot, K,K, K\K,K,K, 
Let us define from (24) that 
n= 7,K,K,K, + t,K, (K, at 137) 
: K,K, (K, ~ T,) (25) 
K,K,(z,+4) 


since K, =1,+7, +. Therefore, 


‘ S 0 if and only if 7, <7). 
2 

Thus, this threshold parameter will be a decreasing function of the 
quarantine parameter T,, resulting in a decrease of disease burden (new 
infections and mortality) if the relative infectiousness of quarantine 
individuals 7, does not exceed its associated threshold 73. By contrast, if 
1), < 1, then the use of quarantine of latent immigrants will increase the 
threshold parameter (R,) and so, increase disease burden. In this case, the 
use of quarantine is detrimental to a community. The result is summarized 
as follows. 


Proposition 1 : The use of quarantine of the latent immigrants will have positive 
(negative) impact in a community if n, S7,. 

Similarly, the rate of change of R, with respect to isolation parameter 
yis given by 


OR, = AIK, (7,K, + a) oF 7,7,K,] 


oy K,K,K,K, 
AIK,K,(7,K, +7,0,)+ 7,K,(K, +737)) (26) 
K,K;K,K, 
We obtain from (26) that 
é K 6, +V, +h 
N 4 a ee) (27) 


O,tV,+u oO, tv, +h 
so that 


OR. 
== 0 if and only if 7, S7,. 
Y 
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It follows that, R, will be a decreasing function of the isolation 
parameter y and hence decrease disease burden if 7, <7,. It will be an 
increasing function of y, on one hand, when 7, > As In such a case, the use 
of isolation will have a detrimental impact in a community. This result is 
summarized in the following proposition. 


Proposition 2 : The use of isolation of infectious individuals will have positive 
(negative) impact in a community if 7, S i 

These results (Propositions 1 and 2) are illustrated numerically in 
section 3.5.1 below. 


100 


Infected Population 


a= 0, 100, 400, 700 


0 50 100-150-200 
Time (days) 


(a) 


an an 
bo 
So 
= 


& 


=) 


b 


Number of Infected humans 
we 

Number of Infected humans 
Se 
S 


S 


0 
200 400 600 800 0 200 400 600 800 
Recruitment rate of humans («) (/day) Recruitment rate of humans () (/day) 


o 


(b) (c) 


Figure 3 


The simulations showing the impact of immigration on the number of infect- 
ed individuals using parameter values as in Table (3). (a) The infected humans 
increases with increasing influx of immigrants. (b) Indicates that the infected 
humans increases to a certain level. (c) Reveals that there is a critical immigra- 
tion rate above which the infected humans decreases. 
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It worthies mentioning here that, the critical parameters 7, and 17; as 
defined in (25) and (27), respectively, remain the same for the case when 
either all immigrants are susceptibles (gq = 1) or the immigration rate is 
negligible, «= 0. Hence, the results in Propositions 1 and 2 also hold true 
for such two cases. 

To explore the impact of influx of immigrants into the susceptible and 
latent compartments on the dynamics of MERS-CoV, we use numerical 
experiment to vary 

different values of the immigration rate (@), while keeping all other 
model parameters fixed. In fact, the numerical simulations reveal that 
the number of infected individuals increases with increasing values of a. 
As depicted in Figure 3(a), when the influx increases from 0 to 700, the 
infected individuals increases from 30 to 90 as well. It also shows that 
the prevalence is very low (ranges between 10 to 14) and the steady state 
value I is also very small, whereas, for 0 < a@< 100, it ranges between 32 to 
42 about, and if a = 700 then 1(200) is less than 100. Moreover, Figure 3(b) 
shows more clearly that the number of infected individuals increases with 
increasing value of a while Figure 3(c) indicates that, there might be a 
critical value of ~@ above which the disease prevalence decreases due to the 
positive impact of quarantine and isolation (Propositions 1 and 2) so that 
the disease will die out with time. However, eradication of MERS-CoV 
in a community is feasible if quarantine and isolation is combined with a 
high reduction of immigration rate. 


Ww 
w 
‘ 


- 
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et 
ae 
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Figure 4 


The simulations showing the impact of quarantine and isolation in controlling 
the infection with parameter values as in Table 3. (a) Quarantine with critical 
value, 7,= 0.3118. (b) Isolation with critical parameter 7, = 0.1087 
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Table 3 
Parameter values for the model (1) 
Parameter Nominal value Reference 

I] 400 per day Assumed 
B 0.25 per day Assumed 
a 20 per day Assumed 
LU 0.0000374 per day [11] 
n, (0.001, 0.65) Assumed 
nN, (0.013, 0.78) Assumed 
1, (0.002, 0.65) Assumed 
Y (0.00138, 0.05) per day [2,8] 

(0.00138, 0.025) per day | [2,8] 
V, (0.0018, 0.0028)per day [2, 25] 
V> (0.0012, 0.0018)per day [2, 8, 25] 
T, (0.0018, 0.0028)per day [2] 
% (0.0012, 0.0018)per day [2, 8] 
6, (0.0337, 0.26) per day [2] 
é, (0.26, 0.5) per day [2] 
p [0, 1] Assumed 
q [0, 1] Assumed. 


3.5.1 Numerical simulations 


The numerical simulations in Figure 4(a) is depicting the result of 
Preposition 1 in which the invasion threshold parameter R,, is a decreasing 
function of the quarantine parameter t,, whenever the modification 
parameter for quarantine is less than a critical value 7, =0.3118. Thus, 
for values of 7, = 0.017, 0.25, less than 73, the cumulative new cases of 
infection decreases from about 27,000 to 16,500 individuals. However, for 
7, = 0.5, 0.75, greater than the critical value, the cumulative new cases 
increases from about 28,500 to 38,500 individuals. 

R, is a decreasing function of the isolation parameter y, whenever 
the modification parameter for quarantine is less than a critical value 
7), = 0.1087. Thus, for values of 7, = 0.001 to 0.017, less than 7, the 
cumulative new cases of infection decreases from about 18,500 to 16,500 
individuals. However, for 7, = 0.25 to 0.55, greater than the critical 
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value, the cumulative new cases increases from about 26,500 to 32,500 
individuals. 


4. Discussion and Conclusion 


A new deterministic model for the spread of the Middle East 
Respiratory Syndrome coronavirus in a community with quarantine and 
isolation as control strategies is formulated and rigourously analyzed. We 
notice that the presented model does not have disease free equilibrium 
point due to constant in flux of the latent immigrants and so, there is an 
endemic persistence of MERS-CoV depending on the invasion threshold 
parameter R, immigration rate as well as the epidemiological status 
of the immigrants. In fact, if the immigration rate a > 0 and fraction of 
susceptible immigrants q € [0, 1), then the disease invades the population 
for all values of R, On the other hand, endemic persistence is only 
possible if either immigration rate is negligible or when all immigrants 
are susceptible provided R; > 1. Stability analysis is carried out for the full 
model and these two immigration scenarios. We showed numerically that 
this equilibrium point is a global attractor when R; <1, see Figure 2. For 
the case when all the immigrants are susceptible, we proved the global 
asymptotic stability of the associated disease free equilibrium (€{) if the 
threshold parameter R; < 1. Such a stability result hold true for the case 
when the immigration rate is negligible, «= 0 and the result follows easily 
from the first scenario. The existence and stability of endemic equilibrium 
for these scenarios follows from that of the full model. 

In order to complete our investigation, we assess the impact of 
quarantine and isolation via a threshold analysis approach. Indeed, we 
showed that quarantine of latent immigrants will decrease the burden of 
the epidemics (new infections), if the rate of the relative infectiousness of 
quarantined individuals does not exceed certain threshold quantity, 17}. 
In a similar note, an isolation of infected individuals using intensive care 
units of hospitals, will have positive impact on the burden of the epidemics 
(new infections and mortality due to disease) by reducing the threshold 
parameter, if the relative infectiousness of isolated individuals does not 
exceed certain threshold quantity, 7;. These results are summarized in 
Propositions 1 and 2 and numerically supported by Figure 4. 

Furthermore, numerical simulations of the model reveal that the 
influx of susceptible and latent immigrants increase the number of 
infected individuals. In fact, one can observe from Figure 3(a) that if the 
influx of these immigrants increases from 0 to 700, the number of infected 
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individuals rises from 30 to 90. Despite such an effect of the immigration 
rate (a), the positive impact of isolation and quarantine drastically reduces 
the disease prevalence. As shown in Figure 3(a), the prevalence is very 
low (ranges between 10 to 14) and the steady state I is also very small. 
While for 0 < a@< 100 it ranges between 32 to 42 about, and if a = 700 then 
1(200) is less than 100. Moreover, Figure 3(b) shows more clearly that the 
number of infected individuals increases as the immigration rate increases, 
whereas, Figure 3(c) indicates that, there might be a critical value of a 
above which the number of infectious humans decreases perhaps, due to 
the positive impact of quarantine and isolation (Propositions 1 and 2) so 
that the disease will die out with time. It is worthy of mention that, most 
of the parameter values used are obtained from some related literature 
as presented in Table 3. The ranges of such parameter values are chosen 
from minimum and maximum as presented in such existing works for the 
sake of consistency. While we assumed the values of the parameters (p, 
q and n, for i = 1, 2,3) in the interval [0, 1] since each of such parameters 
is either a fraction of humans or a fraction that account for the assumed 
reduction in infectiousness of individuals in the latent, quarantined and 
isolated compartments. Moreover, the values of the recruitment rate IT 
and the effective contact rate § are assumed in relation to their values used 
for SARS model in [11]. 

In conclusion, we showed using a deterministic model that MERS- 
CoV can be controlled by quick isolation or monitoring close contacts 
with the patients infected by MERS-CoV and quarantine of suspected 
latently infected individuals. This is similar to the finding in [27]. In 
addition, reducing the influx of immigrants plays a significant role in 
decreasing the number of infected individuals when the recruitment rate 
of immigrants is below a certain critical value. This is a value of a below 
which the prevalence is high. However, community-wide eradication of 
MERS-CoV is feasible when quarantine and isolation are combined with 
high reduction of the immigration rate into a community. 
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